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Abstract 

We propose new techniques for the numerical implementation of the overlap-Dirac 
operator, which exploit the physical properties of the underlying theory to avoid nested 
algorithms. We test these procedures in the two-dimensional Schwinger model and 
the results are very promising. These techniques can be directly applied to QCD 
simulations. We also present a detailed computation of the spectrum and the chiral 
properties of the Schwinger model in the overlap lattice formulation. 



1 Introduction 



The last few years witnessed a major breakthrough in the lattice regularization of Fermi 
fields. The closely related domain wall |TJ and overlap [3| formulations of lattice fermions 
provide a definition of a lattice Dirac operator D which avoids the doubling problem and 
preserves the relevant symmetries of the continuum theory, most notably chiral symmetry 
in the limit of vanishing fermion mass. Unfortunately, this welcome development has come 
at a price: the numerical calculation of the matrix elements of the propagator D~ l and the 
inclusion of Det(D) in the measure entail a substantially increased computational burden, 
which severely constrains the maximum lattice size for viable simulations. Many efforts 
have recently been devoted to finding more efficient ways to perform these computational 
tasks 0. The problem has generally been approached from a numerical analysis point of 
view, looking for approximations that allows one to use sparse matrix techniques, although D 
itself is not sparse. In this article we would like to advocate and explore a different approach, 
where the approximation is based on the physics and proceeds through a projection over a 
subspace, which has a substantially reduced number of degrees freedom but still captures the 
relevant long range properties. Our approximation consists in taking the matrix elements 
of the Wilson Dirac operator in a subspace consisting of the long range modes which we 
expect to dominate the low energy properties of the theory and in constructing the overlap 
Dirac operator numerically, but without further approximations, in this subspace. In the 
complement, the Dirac operator is approximated by its free form. The projection over long 
range modes is done via a Fourier transform after gauge fixing. We also studied a gauge 
invariant method of projection. We will compare the results obtained with our approximation 
in a simplified model with those of an exact calculation, finding satisfactory agreement. It is 
our hope that the approximation may offer a new way to implement the overlap formulation 
in four- dimensional QCD. 

We will focus on the overlap Dirac operator and we will apply our approximation pro- 
cedure to the two-dimensional Schwinger model. We use the overlap formulation because 
it provides an exact framework for the new lattice fermions. The domain wall formulation 
becomes exact in the limit of infinite extent in the extra dimension, in which case it becomes 
equivalent to the overlap formulation up to corrections that vanish with the lattice spacing. 
With finite extent in the extra dimension, the domain wall formulation can be viewed as 
a numerical approximation to the exact result. Since one of our goals is to compare the 
results of our approximation to exact results, the overlap formulation is better indicated. 
We study the Schwinger model because the system is simple enough that we will be able to 
calculate D~ l and Det(-D) exactly as well as in our proposed scheme of approximation. The 
results of the exact calculations are, we believe, interesting per se, because they are more 
extensive than what, to the best of our knowledge, has been done up to now and validate in 
an impressive manner the advantages of the new formulation of lattice fermions. 
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The plan of this paper is as follows. In the next section we will briefly review known 
properties of the Schwinger model in the continuum and in the lattice overlap formulation. 
In Section |3] we will present the results of a numerical simulation where propagator and 
determinant of the overlap Dirac operator are calculated exactly. In Section [| we will intro- 
duce the proposed scheme of approximation, which is based on the projection over a reduced 
number of degrees of freedom. We will compare the results obtained with this approximation 
with those of the exact calculation. In Section |S| we will describe an approximation to the 
determinant of D, based on a coarsening of the lattice and compare results obtained with 
the full and approximated determinants. In the last section we will present a few words of 
conclusion. 



2 The Schwinger Model and the Overlap Formulation 

We consider the Schwinger Model in (two-dimensional) Euclidean space-time and we allow 
for Nf degenerate flavors. The model is defined by the action 



S = d 2 x 



A F[iyF[xv 

i=l,N f 



(1) 



where ipi, ipi are two components spinors, is the U(l) gauge potential, = d^ + igA^ and 
F^v = d^Ay — d v A^. We will use the following two-dimensional representation of 7-matrices 

7i = °1 ) 72 = 02 , 75 = -«7i72 = °3 (2) 
where Oi are the Pauli matrices. 

The topological charge of a classical gauge configurations is 

Q(A) = £-jd 2 xe^F fMU (x) (3) 

while the index of the Dirac operator is given by the difference between the numbers of 
positive (n + ) and negative (n_) chirality zero modes 

index(y4) = n_ — n + (4) 
The Atiyah-Singer theorem states that 

Q(A) = index(v4) = n_ - n + (5) 
Moreover in two dimensions the vanishing theorem ensures that 

n + 7^ — > n_ = 

n_ ^ -> n + = (6) 
2 



At the quantum level the Schwinger model does not require infinite renormalization: g and 
m are finite bare parameters. In the massless limit the model can can be solved exactly f|. 
In the following we will analyze the systems with Nf < 2. Information for generic Nf can 
be found in 0. 



2.1 The Schwinger model with Nf = 1. 

The classical theory has a vector U(l)v symmetry and a softly broken axial symmetry U(1)a, 
which at the quantum level is also broken by the anomaly The Ward Identities (WI) for 
the corresponding quantum theory are 

d„V„ = , d„A„ = ^-e^ v + 2mP (7) 

where the Axial and Vector currents are defined as 

V n( x ) = ^( x h^( x ) > A ^( x ) = 4>( x hvl5^( x ) = -ie^K , (8) 
and the scalar and pseudoscalar densities are 

S(x) = , 4>{x)'4){x) , P(x) = ■0(x)75^(a;) . (9) 

It is a peculiar property of the two-dimensional space-time that the vector and the axial 
currents are not independent of each other. 

In the massless case a free fermion field factorizes and the Schwinger photon acquires a 
mass due to the anomaly. The photon's field $ is 

V» = -^e, u d u $ (10) 

'7T 



and the vector current correlator 

(V,(x)V u (y))o (ll) 

is the same as for a free massive propagator with mass (/^o/g) 2 = 1/vr. The chiral condensate 
is 

g ~ 2tt3/2 

where 7 = 0.577216 ... is Euler's constant. Note that the formation of the fermion conden- 
sate, as in one-flavor QCD, does not imply spontaneous symmetry breaking; the £7(1) chiral 
symmetry is already broken by the anomaly. 

For small masses the corrections to the massless results can be obtained using chiral 
perturbation theory || |9j and the Schwinger mass at the first order is given by 

(m 2 ,l_ 47r M£p) (is) 
v 9 ' n 9 v 9 J 



2.2 The Schwinger model with Nf = 2. 



In the model with two degenerate massive flavors, the classical theory has a U(l)y x SU(2)y 
vector symmetry, an axial SU(2)a symmetry softly broken by the mass term and a U(1)a 
symmetry broken by the quantum corrections. The Ward identities corresponding to the 
axial and vector £7(1) symmetries are analogous to the previous case. The Ward identities 
associated to the non-singlet axial and vector currents are given by 



W = (14) 



and 
where 



d^A; = 2mP a (15) 

^ a = ^y7^, ^ = ^y7^, P a = ^y7s^ (16) 

and A a = a a are the generators of the SU(2) in flavor space. The non-singlet axial Ward 
identity ( |T5| ) will be one of the key ingredients to test the properties of the overlap regular- 
ization. 

In the massless limit, there are one massive singlet particle (r/) with mass (/j, v /g) 2 = Nf/n 
and a triplet of massless particles (it). Analogously to the Nf = 1 model, the correlation 
functions of the vector currents are the same as for free particles. Unlike QCD, the chiral 
condensate (ipip) o = 0. In fact a non-zero fermion condensate would break spontaneously the 
SUa{2) chiral symmetry of the model. But spontaneous symmetry breaking is not possible 
in two dimensions HTUI. In this case also, the corrections to the massless results can be 



obtained by a semi-classical analysis ]TTJ and the masses for triplet and singlet areQ 

2 5/6 





( m 
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- -+( 


9 J 


71 ' 




e 2 ^ 3 —- = 2.163... (17) 



2.3 The Dirac Operator in the Overlap Formulation. 

In the lattice regularization of the Schwinger model, the fermionic fields are defined over the 
sites of a square lattice, with lattice spacing a, and the gauge potentials are replaced with 

1 Recently a more precise computation in the limit of large coupling g and small mass m has been performed 
in Q. It gives c = 2.008 . . . 
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£7(1) link variables U^x) defined over the oriented links of the lattice. The Euclidean lattice 
action is given by 



N, 



Sl=I3 £ [l-ReU^(x)]+a 2 J2J2^( x ) 

x,[i<v i=l x,y 



, am 

1 - —)D(x,y) +m 



My) 



where (3 = l/(ag) 2 , g being the bare coupling constant, 

U^(x) = U^(x)U v (x + afl)Ul(x + au)Ul(x) 



(19) 



(x being a unit vector in direction /i, and D is the lattice Dirac operator in the overlap formula- 
tion, as introduced by Neuberger. Occasionally we will also refer to D as the Neuberger-Dirac 
operator. D is constructed as follows. One starts from the massless Wilson-Dirac operator 



M M 



(20) 



where and V* are the forward and backward lattice derivative, i.e. 



V^{x) 
V>0r) 



~[Un(x)if)(x + afi) - i>{x)} 
I r 

- ift(x) — UUx — ap,)ij}(x — a/2) 



(21) 



One performs a polar decomposition of D w — 1/a, expressing this operator in terms of a 
unitary operator V and its modulus 



D w = V 

a 



Dl - -) (D w - - 

aj \ a 



(22) 



From Eq. it follows that V is given by 



V 



D 



w 



D 



u ) (Dw 

ay \ a 



Finally, the Neuberger-Dirac operator for a fermion of bare mass m is given by 

1 m s 



D 



[l + V)+m 



(23) 



(24) 



The Neuberger-Dirac operator satisfies the 75-Hermiticity condition 

£>t = 75^75 



(25) 
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Moreover, for m = it satisfies the Ginsparg- Wilson (GW) relation 
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75/T 1 + zrSs = a 75 (26) 
which implies that the fermion action (|T8"D at finite lattice spacing is invariant under the 



following continuous symmetry [14 



<ty = 7b(1 - a£)V>, 64> = $y 6 (27) 

This can be interpreted as the lattice form of chiral symmetry. The corresponding flavor 
non-singlet chiral transformations are defined including a flavor group generator in Eq. ( |2"7| ) . 

The 75-Hermiticity condition fl25|) and the GW relation ( p6|) imply further algebraic 
relations which can turn out to be useful in numerical simulations Ul5| . In particular 



D + D ] = aDD ] = aD^D (28) 

hence D and commute, i.e. D is normal and therefore the eigenvector are mutually 
orthonormal. 

It is important to stress that locality, the absence of doubler modes and the correct 
classical continuum limit are not guaranteed by the GWR in Eq. fl2"Bj). Indeed, there exist 
lattice fermion actions which satisfy the GWR but which do not meet the above requirements 
]16| | . The Neuberger operator satisfies all the above requirements and is local in the weak 



coupling regime, as shown in Ref. |17] . 

The geometrical definition of the topological charge on the lattice is 

Q = -Llm$>g(C/ 12 ) (29) 



The flavor- singlet chiral transformations in Eq. fl27|) leads to chiral Ward identities analogous 
to the ones in Eq. (0) and the anomaly term arises from the non-invariance of the fermion 
integral measure [|14| . The non-singlet chiral transformations analogous to Eq. (|27|) and the 



locality of the Neuberger-Dirac operator imply axial Ward identities analogous to Eq. ( |T5D 
and therefore the quark mass renormalizes only multiplicatively, i.e. the critical bare mass is 
zero, the corresponding conserved axial and vector currents do not need renormalization and 



there is no mixing between 4-fermion operators in different chiral representations |T^, [19 
Finally, the chiral condensate requires a subtraction and is defined as 

X = ~m - \D)^) (30) 
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3 Numerical Results from an Exact Calculation 



In order to test our method of approximation and, at the same time, to increase the body 
of information on the lattice Schwinger model, we performed an extensive simulation with 
the exact calculation of Neuberger-Dirac operator and its inverse. For previous numerical 
work on the Schwinger model in the overlap regularization see [P0|-||24||. We considered 



the Schwinger model with (3 = 6 on a square lattice with N x = N y = 24. On a lattice 
of this size, the discretized Dirac operator is a complex matrix of dimension 1152 x 1152, 
for which we could use full matrix algebra subroutines without excessive burden on the 
resources available to us (the whole calculation used approximately 20,000 processor hours 
on the Boston University SGI/Cray Origin 2000 supercomputer). Having settled on the 
lattice size, we selected (3 = 6 on the basis of previous results which indicated that, for most 
of the fermion masses we were planning to consider, the lattice would span at least a few 
correlation lengths. We generated 500 configurations of the gauge variables U^{x) distributed 
according to the pure gauge measure 

Sg = P E [1-RetWx)] ( 31 ) 

X,fl<U 

(see also Eq. [T9| ) . These were obtained by a standard Metropolis Monte Carlo simulation, 
with 10000 upgrades of the whole lattice between subsequent configurations. With the 
pure gauge measure the plaquettes are essentially uncorrelated, apart from the constraint 
coming from the periodic boundary conditions, so the procedure we followed should be amply 
adequate to produce independent configurations. For calculations with one and two flavors 
of dynamical fermions, we incorporated the determinant of the lattice Dirac operator in the 
averages giving the observables. While with large variations of the determinant this way of 
proceeding would lead to an unacceptable variance, in the present calculation we found the 
range of values taken by Det(D) to be sufficiently limited to warrant our averaging procedure 
(see Fig. [L^ in Section [5].) This is of course due to the rather small size of our system. With 



larger systems one should incorporate the determinant (or a suitable approximation to it) 
in the measure used for the simulation. 

For each configuration, we performed a singular value decomposition 

D w -- = UAU (32) 
a 

where U, U are unitary matrices and A is diagonal, real and non-negative. The operator V 
in Eq. ^ is given by 

V =UU (33) 

We proceeded then to the diagonalization of V, calculating all its eigenvalues and eigen- 
vectors. From these, it is clearly straightforward to calculate both the determinant of the 
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Figure 1: M w /g vs. (m/g) 2 / 3 for the full operator for Nj = 2. The dashed line represents 
the fit as in Eq. of the four highest masses. 



Neuberger-Dirac operator D as well as its associated propagator D^ 1 for any value of the 
fermion mass m. The full diagonalization of V is computationally more demanding then 
the direct calculation of D~ l , which typically gives also Det(D) as a by-product, but we 
were interested in the actual spectrum of V for comparison with the approximations that 
will be discussed later. From the fermion propagators we calculated the meson propagators 
projected over zero momentum. We focused on the vector correlators 



( J2 ^( x > yh2ip(x, y)*P(x + t, y'h 2 ip(x + t, y')) 



(34) 



because they are saturated by a single particle contribution in the massless limit. Practi- 
tioners of lattice calculations will certainly appreciate the value of being able to sum over 
all source locations, as opposed to having to deal with selected columns of the meson prop- 
agators. Of course we added to the averages the correlators obtained from the interchange 
of x and y in Eq. |33| for a further gain in statistics. 

From fits to the meson correlators we extracted the meson masses in a standard manner. 
Figures [I], 0, 0, f| reproduce our results for the meson masses as functions of the fermion 
mass. The values we obtained for the masses are also reported in the tables included in 
Sect. The errors have been evaluated with the jackknife method. Figure [I] illustrates the 
behavior of the isotriplet mass M n in the model with two flavors. This mass is expected 
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Figure 2: m AWI / g vs. m/g for the full operator for Nj = 2. 

to vanish for m = and chiral perturbation theory predicts an m 2 / 3 dependence on m (see 
Eq. [17p. The dashed line in the figure shows the results for the fit 

Mjg = A + B(m/g) 2 / 3 (35) 

for the four heaviest masses, which satisfy the condition N X M V > 4 and therefore are less 
likely to be affected by finite volume effects. This fit gives A = —0.001(65) and B = 2.10(14). 
We also performed the fit 

Mjg = C(m/ g y (36) 

which gives C = 2.10(17) and 7 = 0.67(7). These results confirm in an impressive manner the 
chiral properties of the Neuberger operator and the mass dependence expected by analytical 
calculations O . 

Equally gratifying is the comparison of the value of the fermion mass m in the Lagrangian 
with the value m AWI that can be extracted, up to discretization effects, from the lattice 
analog of the axial Ward identity in Eq. |15[ As shown in Fig. |2|, the intercept and the 
slope of the linear fit through the four heaviest masses are compatible with zero and one 
respectively. 

In Fig. |3] we illustrate the behavior of the singlet mass M n versus m 4 / 3 , always for Nf = 2. 
The singlet meson propagator is given by the difference of the connected and disconnected 
terms in the correlator and, because of the cancellations, the errors are larger than in the 
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Figure 3: M v /g vs. (m/g) 4 / 3 for the full operator for Nf = 2. 



triplet case. The numerical results are, however, consistent with the theoretical prediction 
of Eq. |T3. 

Finally, in Fig. |] we display the singlet mass in the single flavor case. The errors are 
smaller than in the two flavor case because of the different relative weight of the connected 
and disconnected contributions. Again we find reasonable agreement between numerical 
results and the theoretical prediction. 



4 Approximation of the Neuberger-Dirac Operator 

Anybody who has ever seen a plot of the spectrum of eigenvalues of the Wilson-Dirac operator 
cannot but be left with the impression that there is a huge redundancy of states. The 
projection over a translated unitary circle done by the Neuberger-Dirac operator avoids the 
problem of mode doubling, but preserves the overall count of states. Yet, it would appear 
that the physical properties of the system should be determined by the eigenstates of D with 
eigenvalues in the neighborhood of A = m (i.e. the eigenstates of V with eigenvalue Xy closest 
to —1, cfr. Eqs. ^3], p4|) , since these are the states with a smooth long range behavior, expected 
to go over the physical states of the continuum in the limit a — > 0. This suggests that it 
should be possible to reconstruct the physical observables from these states alone. The way 
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Figure 4: M v /g vs. (m/g) for the full operator for Nf = 1. 



in which such states contribute to physical observables has been studied in the literature ||21 



Here we would like to make the point that the overlap formulation is particularly well suited 
for the implementation of the above approximation, since the unitarity of V provides, in some 
sense, a checkpoint on the approximation itself. Of course, one must be careful in attempting 
any approximation based on neglecting the short-wavelength part of the spectrum, even if 
the corresponding states are largely lattice artifacts, since one knows that in quantum field 
theory the infrared and ultraviolet components of the spectrum are subtly related. Thus, the 
chiral eigenstates with Ay = —1 ("zero modes"), which V exhibits in presence of a gauge field 
with non-trivial topology, find their counterpart in states with opposite chirality at Ay = 1. 
Here too, however, the special features of the overlap formulation come to the rescue, since, 
as we will show, both the presence of zero modes and the correspondence between Ay = — 1 
and Ay = 1 eigenstates of opposite chirality can be preserved by the projection over a subset 
of physical states. In this section we will illustrate such a projection by comparing the 
approximate results for the observables with those of the exact calculation. 

One possible scheme of approximation which is computationally very convenient consists 
in performing a projection over states of low momentum in Fourier space, after gauge fixing 
to a smooth gauge field configuration. In a smooth gauge, because of the suppression of 
short-wavelength fluctuations due to asymptotic freedom, one expects the structure of the 
Wilson operator come more and more diagonal in momentum space for increasing momenta. 
This notion is schematically illustrated in Fig. [5] and also underlies the technique of Fourier 
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acceleration for the calculation of quark propagators [25]. Accordingly, we implemented the 
following approximation. 

D(q,p) 





Figure 5: Structure of the Wilson operator in momentum space and in a smooth gauge. 



We fixed the gauge on all configurations to the Landau gauge by demanding that the 
function 

G = J2MU,(x)] (37) 

be maximal. A relaxation procedure produces several local maxima (Gribov copies). Among 
all these configurations, we selected the one that produced the absolute maximum for G, 
subject to the further constraint that for all x and li ReU^x) > 0.5. Some care must 
be exerted with the configurations that have non- vanishing topological number Q. These 
configurations cannot be brought to a uniformly smooth gauge. To discuss them further it 
is convenient to introduce the notation 

U^x) = e l9 ^ x) (38) 
U^(x) = e^»W (39) 

Then, with Q ^ there must be one or more plaquettes where 

dfj,u( x ) — + ®v{ x + fl p0 — 0fi{x + au) — 9 v {x) + 2mn (40) 
with m a non- zero integer. This implies that for such plaquettes 

Abs[^(x) + 8 v (x + a fx) - O^x + au) - 9 v {x)\ > n (41) 
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and thus the corresponding gauge variables cannot be all close to unity. We will call these 
plaquettes the "hot spots" of the gauge configuration. Their location is gauge dependent 
but, as long as Q 7^ 0, they cannot be eliminated. For the configurations Cq with Q ^ 
we fixed the gauge as follows. First we created a configuration C with the same value of Q 
by superimposing gauge field configurations in the Landau gauge with uniform U^ u {x) and 
Q = 1. We placed the hot spots of this configurations in the locations where the action of 
the original field configuration was smaller, using a weighting procedure that minimized the 
overlap of hot spots. Then we brought to the Landau gauge the gauge field configuration 
C" obtained subtracting C from Cq. Finally we superimposed C to the gauge transformed 
C" . The resulting configuration is gauge equivalent to the original Cq, is approximatly 
in the Landau gauge and is "reasonably smooth", the latter statement being justified a 
posteriori by the application of our approximation. Clearly the procedure described above 
relies heavily on the Abelian nature of the gauge field and will have to be generalized to 
handle non-Abelian systems. 

After gauge fixing, we project the lattice Wilson-Dirac operator over the subspace T of 
Fourier space spanned by the eigenvectors with lowest momenta. For the results we present 
in this paper, the projection was done over the subspace defined by 

M 2 + {p y af < {p c af (42) 

where the cut-off momentum p c was chosen so that the total number of momenta included 
in the projection equals 145, i.e. a fraction z ~ 1/4 of the total number of momenta. (We 
could not choose p c in such a way as to project over exactly one fourth of the space, since 
the number of momenta satisfying Eq. is always odd.) Accordingly, the dimensionality 
of T is 290, out of a total dimensionality of the fermion field equal to 1152 (24 x 24 lattice 
sites x2 spin degrees of freedom). This corresponds to an effective reduction of a factor 
of two for each dimension of the lattice. As a matter of fact, we experimented with even 
smaller values of z and found that, with Q = 0, z can be chosen substantially smaller than 
1/4 without spoiling the good quality of the approximation. For configurations with a non- 
trivial topology the infrared properties of the Wilson operator are still well reproduced on 
the subspace. However the real doubler modes of the truncated Wilson operator are shifted 
towards the infrared region and, as they pass the projection point, the association between 
topology and chirality is lost. 

The approximation to the Neuberger-Dirac operator is obtained following the construc- 
tion of Eqs. ^ and where D\y and V are replaced by their projections over T: D\y and 
V (see also |2£j). If we denote by D$ the projection of the free Dirac operator over the 
complement of T and by V"*- ) the corresponding unitary factor in its polar expansion, we 
finally take 

V approx = V® V® (43) 
as our approximation to V. More specifically, given any fermionic vector ip, we project it over 
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T and its complement: tp = ipp + ^F- Vapproxtp is then given by Vipp + V^ .The approximate 
form of the Neuberger-Dirac operator follows from Eq. 




Figure 6: Comparison of spectra on a configuration with Q = 1. 



The first check on the approximation is that it should reproduce well the spectrum of 
long range modes, i.e. those with V eigenvalues in the neighborhood of -1. In particular, 
configurations with Q ^ should have \Q\ eigenvectors with V eigenvalue exactly equal to 
-1. This turned out to be the case. With the chosen cut-off p c (See Eq. [|2p the approximation 
consistently gives satisfactory results for the spectrum of long range modes. We display in 
Fig. the spectrum of eigenvalues of the exact 1 + V (x symbols - blue) and the approximate 
1 + V (plus symbols - red) for a configuration with Q — 1. The matching in the physical part 
of the spectrum, close to 0, is excellent (see the blow-up at the right of the figure). Also, 
as we anticipated above, the approximation preserves the presence of a zero mode, which 
is a chiral eigenstate. It is worthwhile to notice that V is constructed from a projected 
Wilson operator D\y which has the same 75 transformation properties as the original Wilson 
operator. Thus V will have the same 75 transformation properties as V and, in particular, its 
eigenvectors with eigenvalue —1 will be matched by a corresponding number of eigenvectors 
with eigenvalue 1. Since V® has no eigenvectors with eigenvalues ±1, this carries over to 

the V a pp rox - 

Closely related to the spectrum is the value of the chiral condensate. It is instructive to 
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Figure 7: Chiral condensate from the exact calculation. 
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Figure 8: Chiral condensate with the Fourier approximation. 
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reexpress the condensate of Eq. |30] in terms of the eigenvalues of V. With a = 1 we find 



x = -j r $Q--\D)4>) 



(Det( J D) Af /Tr 



D- 



N x N y l-m/2 



(Det(D) 



A' 



(44) 



(We use the subscript U in the r.h.s. of this equation as a reminder that the averages are 
quenched averages over the gauge field. See also the discussion following Eq. [33].) Using 
Eq. |4| we get 



1 (Det(D) N f ^ + fJ2 



X 



1— cos < 



'3 (l+rn 2 /4)+(l-rn 2 /4)cos(f> J 



U 



N X Ny 



(Det(D) N f) 



u 



(45) 



where we separated in the sum the \Q\ pairs of eigenvalues equal to -1 and 1 from the others, 
which occur in pairs of complex conjugate values exp(±z^). 

Figure ^ reproduces the result of the exact calculation, Fig. || the result obtained by 
inserting in Eq. f|5] the eigenvalues obtained with the Fourier approximation. The two sets 
of results are indistinguishable (we did not plot them on the same graph, because one set of 
lines would have covered the other). The crosses on the two figures show the analytic result 
for to = 0, Nf = 1. 
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Figure 9: Results for ^y 2 -triplet propagator on a configuration with Q = 0. 



Another test of the approximation is obtained by comparing configuration by configura- 
tion the propagators for fermion bilinears. Our approximation gives quite satisfactory results 
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Figure 10: Results for ^-triplet propagator on a configuration with Q = — 1. 



20 



16 



12 



U 



4 - 



Full Op. 

FT Op. I 



4 



4 



4 



12 
t 



If 



24 



Figure 11: Results for ^-triplet propagator on a configuration with Q = 4. 
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Figure 12: Disconnected part of the y 2 -singlet propagator on a configuration with Q = 0. 
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Figure 13: Disconnected part of the y^-singlet propagator on a configuration with Q = — 1. 
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Figure 14: M^/g vs. (m/g) 2 / 3 for the full operator and the gauge-fixing approximation for 
Nf = 2. 



for the connected Green's functions. We reproduce in Figs. |9|, |1(] and |11] the propagator of 
the 72-triplet bilinear in three randomly chosen configurations with Q — 0, — 1 and 4, respec- 
tively. We see that propagators obtained from the approximation to the Neuberger-Dirac 
operator compare quite well with those obtained from the exact operator, even with Q = 4. 

Our results for the disconnected Green's functions are not conclusive. We reproduce in 



Figs. 12, 13 the disconnected parts of the propagator in the same configurations as used 



for Figs. H], 10. The approximate propagators are smoother than the exact ones, which 
can be understood as an effect of the truncation over the short wavelength modes, but 
not incompatible in an average sense. Of more concern is that we found the fluctuations 
of the approximate disconnected propagators to be substantially larger than in the exact 
calculation. The large fluctuations of the disconnected components in the exact calculation 
already makes the error in the singlet masses rather big. The even larger fluctuations of the 
disconnected components in the approximate calculation prevented us from obtaining in this 
case a meaningful result for the singlet masses. 

We reproduce in the tables all the masses which we were able to calculate. We see that 
whenever we could get a mass value from the Fourier approximation, the result turned out 
in reasonably good agreement with the exact calculation. A fit to Eqs. |35] and |36] of the 
approximate calculation gives the following results for the parameters: A = —0.005(65), 



19 



B = 2.10(15), C = 2.10(17) and 7 = 0.67(7) with excellent agreement with the results of the 
full calculation, presented immediately after Eqs. [35] and [36[ A comparison between masses 
in the exact and approximate calculations is also presented in Fig. [H. 



P = 6.0, V = 24 2 


m/g 


m n /g mf wl /g 


Full Operator 





0.15(1) -0.0026(2) 


0.0244 
0.0485 
0.0960 
0.1427 
0.1884 
0.2333 


0.393(13) 0.0234(1) 
0.442(10) 0.0472(2) 
0.553(8) 0.0949(3) 
0.661(6) 0.1420(4) 
0.761(6) 0.1884(4) 
0.857(5) 0.2343(5) 


Fourier Approximation 


0.0 


0.14(1) -0.0049(3) 


0.0244 
0.0485 
0.0960 
0.1427 
0.1884 
0.2333 


0.385(13) 0.0208(3) 
0.435(10) 0.0427(5) 
0.548(8) 0.0875(8) 
0.656(7) 0.1321(11) 
0.758(6) 0.1761(13) 
0.855(5) 0.2194(16) 



Table V. Masses for Nj = 0. 



In closing this section let us mention a different, gauge invariant scheme of approximation, 
which we also tried. We built the subspace T used for the projection of Eq. ^3 out of the 
lowest eigenvectors of the negative of the covariant Laplacian, namely the operator (see 
Eq.© ' 

-A = -£V M V; (46) 

The lowest eigenmodes of —A are typically slowly varying, whereas the highest modes ex- 
hibit rapid variations from site to site. Thus the subspace JF formed out of the lowest 
eigenvectors of —A is well suited to represent the physical excitations of the system. The 
construction of this subspace is gauge-invariant and, contrary to the approximation based 
on the Fourier transform, does not require any gauge fixing. Finding a sufficiently large 
number of eigenvectors of —A is however computationally expensive. An even more serious 
problem is that this approximation does not provide any simple way to define the projection 
of a "free operator" on the complement of T: the form of the free Wilson operator is clearly 
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(3 = 6.0 


, V = 24 2 


m/g 


m n jg 


Analytic 





0.5642 


Full Operator 





0.67(6) 


r\ no A A 

0.0244 


0.74(10) 


0.0485 


0.77(8) 


0.0960 


0.83(5) 


0.1427 


0.90(4) 


0.1884 


0.98(3) 


0.2333 


1.05(3) 


Multig 


;rid Det. 





0.68(6) 


0.0244 


0.75(10) 


0.0485 


0.78(7) 


0.0960 


0.84(5) 


0.1427 


0.91(4) 


0.1884 


0.98(3) 


0.2333 


1.06(3) 



Table 2: Masses for N* = l. 



gauge dependent. Thus, for lack of better alternatives, we replaced V(0) in Eq. |3] with the 
identity operator, projecting all the remaining eigenvalues to 1. As a consequence of this 
rather drastic projection, the propagators of the bilinears turned out to exhibit oscillations 
around the exact propagators, with much less satisfactory agreement than we found in the 
Fourier approximation. While these oscillations do not necessarily rule out the possibility 
of still obtaining good values for the masses through a suitable fitting procedure, the worst 
quality of the bilinear propagators together with the higher computational costs dissuaded 
us from pursuing this approximation further. 

5 Coarse grid Approximation to the Determinant 

In the previous section we explored the possibility to truncate the Neuberger operator to its 
long range components on a given gauge configuration. In the spirit of long range approx- 
imations, there is another possibility that can be explored, which is the coarse graining of 
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f3 = 6.0, 


V = 24 2 




m/g 


m AWl 


m n /g 


m v /g 


Analytic 











0.7979 


Full Operator 


0.00 


-0.004(2) 


-0.001(65) 


1.00(25) 


0.0244 


0.0233(2) 


0.14(3) 


1.0(5) 


0.0485 


0.0468(6) 


0.26(5) 


1.0(4) 


0.0960 


0.0945(13) 


0.44(4) 


1.1(2) 


0.1427 


0.1420(14) 


0.57(3) 


1.1(2) 


0.1884 


0.1890(13) 


0.69(2) 


1.2(1) 


0.2333 


0.2353(11) 


0.79(1) 


1.23(9) 


Fourier Approximation 


0.00 


-0.002(2) 


-0.005(66) 




0.0244 


0.018(1) 


0.11(3) 




0.0485 


0.044(1) 


0.24(5) 




0.0960 


0.092(1) 


0.43(4) 




0.1427 


0.139(1) 


0.57(3) 




0.1884 


0.184(1) 


0.68(2) 




0.2333 


0.228(1) 


0.79(1) 




Multigrid Determinant 


0.00 


-0.004(1) 


-0.002(40) 


1.09(25) 


0.0244 


0.0234(1) 


0.14(2) 


1.1(4) 


0.0485 


0.0469(4) 


0.26(3) 


1.1(4) 


0.0960 


0.0946(9) 


0.44(2) 


1.1(2) 


0.1427 


0.1420(10) 


0.57(2) 


1.2(2) 


0.1884 


0.1889(9) 


0.68(1) 


1.2(1) 


0.2333 


0.2351(8) 


0.79(1) 


1.26(9) 



Table 3: Masses for Nt = 2. 



the Dirac operator, as done in multigrid calculations [27 



A preliminary step is the projection of the gauge field to a coarser lattice. Starting 
from an arbitrary gauge configuration on an N x N lattice, we apply a blocking procedure 
to arrive at a gauge configuration on a coarser A^/2 x N/2 lattice designed in such a way 
that the coarser lattice should carry the long range features of the finer N x N lattice. The 
absence of additive mass renormalization then allows one to define a Neuberger operator on 
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Figure 15: Sketch of the blocking procedure, 
coarse lattice are fixed to local Landau gauge. 
Eq. £^ to give links of the blocked lattice. 



Original Links corresponding to a site in the 
Links between them are averaged according to 



the coarser lattice, which is related to the Neuberger operator on the finer lattice without 
tuning of the bare mass and reproduces all its long range features. One could then in 
principle measure the Green's functions on the coarser lattice. It would however be no 
trivial problem to interpolate them back to the finer lattice compensating for the errors 
introduced by coarsening. We therefore decided only to measure the determinant on the 
coarse lattice and compare it to the one on the fine lattice. 

To block a gauge field U^(x, t) on an N x N lattice down to a gauge field U^x, t) on an 
N/2 x N/2 lattice we applied the following procedure: 

• Divide the sites of the fine lattice into N/2 x N/2 blocks of 2 2 sites. The center of each 
of these blocks will correspond to a site (x, t) in the coarse lattice. 

• Fix the gauge such that the links connecting sites within one block are as close as 
possible to unity, i.e., 

JV/2-1 

{U & {2x, 2t) + U { (2x, 2t) + U £ (2x, 2i + 1) + Uf(2x + 1, 2t)) max (47) 

x,t=0 
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Figure 16: logDet(D) on the blocked lattice versus logDet(-D) on the full lattice, for lightest 
mass, m/g — 0.0244, on the full system. 
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Figure 17: M n /g vs. (m/g) 2 / 3 for the full determinant and for the blocked one. 
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Construct the link between two sites of the coarse lattice as close as possible to the 
links between the corresponding fine lattice blocks, i.e., 



((Uj(2x, 2t + 1) + U\(2x + 1,2*+ l))Ui(x, t) 
((Ul(2x + 1, 2t) + Ut(2x + 1, 2i+ l))f4(x, t) 



max and 



max 



(4i 



Figure [L5| gives a picture of how the blocking is done. 

As stated above, we calculated the fermion determinant on the blocked lattice and com- 
pared it to the determinant on the full lattice. As shown in Fig. [16|, the blocked determinant 
follows the full determinant closely. Moreover the eigenvalues of V exactly equal to —1 are 
always preserved by the blocking. Thus one would expect the effects of dynamical fermions 
to be well reproduced by an approximation where the blocked determinant is used instead 
of the full determinant. 
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Figure 18: M v /g vs. (m/g) 4 / 3 with the full, blocked and quenched determinant for the full 
operator for Nf = 2. 



This is confirmed by results shown in Figs. [17] and |18| In Fig. [17| we reproduce the 
values of the isotriplet mass M n obtained by using either the full determinant or the blocked 
determinant in the calculation of the observables. We see that the results obtained with the 
blocked determinant are in very good agreement with those obtained with the full determi- 
nant. The physical results are those for Nf = 2. We inserted in the figure also the results 
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of a quenched calculation and of a calculation with Nf = 1 in order to highlight the effects 
of the determinant. In particular, one sees that the data support a linear dependence of 
M w versus (m q /g) 2 ^ 3 with M w vanishing in the chiral limit only when one has 2 dynamical 
quark flavors, i.e. when one uses the square of the full determinant Det 2 (.D) as weight factor. 
(This observation is based on the data for the 4 rightmost points in the figure, which have 
N x M n < 4 and thus are less likely to be affected by finite size effects - see also the fit in 
Fig. [T]and the discussion that follows.) This nontrivial result does not change however, when 
Det 2 (.D) is replaced by the square of the blocked determinant Det 2 lock (D). 

We reproduce in Fig. [18] the results for the singlet mass m v for Nf = 2. Once again we 
also include the results of a quenched calculation to show the effect of the determinant. In 
this case also unquenching with the blocked determinant is consistent with unquenching with 
the unblocked determinant. The results however are less impressive, since the observable 
contains disconnected propagators and is therefore very noisy. 

6 Conclusions 

We have explored, in the context of the overlap formulation of lattice fermions, a scheme of 
approximation which focuses on the physical properties of the system and takes advantage of 
the special features of the Neuberger-Dirac operator. Our method of approximation is based 
on the projection of the Wilson operator on a subspace of lower dimensionality, where the 
Neuberger-Dirac operator and its inverse are constructed by dense matrix techniques. We 
tested our approximation on the Schwinger model, finding satisfactory results. Of course, 
the main interest of any scheme of approximation for overlap fermions is in its applicability 
to four- dimensional QCD. The possibility of a successful extension of our method to this 
theory will depend on how much one can reduce the dimensionality of the space used for the 
projection with respect to the full space, while keeping a good approximation to the physical 
observables. 

In the two-dimensional Schwinger model, the simplest implementation of these ideas 
has been able to reduce the dimensionality of the fermionic vector space by a factor of 
approximately 4. If this is an indication that one can achieve a reduction of a factor of two 
per space-time dimension, it would mean that in four dimensions the dimensionality of the 
fermionic space can be reduced by a factor of order 16. This may not be enough for QCD 
calculations on lattices of realistic size, but even larger reductions of dimensionality are not 
necessarily out of the question. In our experiment with the Schwinger model, we found 
that we could reduce the dimensionality of the subspace used for the projection by a factor 
ranging well above 4 (as much as 10 or more while still keeping a good approximation to 
physical observables) for the configurations with trivial topology. As mentioned in Section |], 
for configurations with non-trivial topology the infrared properties of the Wilson operator are 
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well captured even on small subspaces, but the Neuberger projection mixes then eigenvectors 
with different chiralities. This could however be more a shortcoming of the basis used for 
the projection, i.e. the Fourier basis, than of the method of approximation in itself. 

Even if it turned out that the reduction of dimensionality one can achieve is not sufficient 
for a practical use of dense matrix techniques in the projection subspace, the approximation 
we have studied can still be of value. The Neuberger-Dirac operator and its inverse could 
be calculated within the subspace with sparse matrix numerical techniques similar to those 
currently in use for the whole space, with the advantage of having to deal with a substantially 
smaller system. Alternatively, the approximate eigenvectors obtained by the projection could 
be used for a preconditioning of the calculation of the propagators on the full space. This 
is a possibility which we have not explored in our work, mostly for limitations of time and 
resources, but which might be worth studying in a future investigation. It is likely, though, 
that incorporating in the computational techniques used for dealing with overlap fermions 
as much insight as possible on the properties of the physical excitations will pay handsome 
rewards, and we are currently investigating the application of our approximation method to 
four-dimensional lattice QCD. 



Acknowledgments 

We wish to thank S. Capitani, H. Neuberger and U. Heller for interesting discussions. This 
work has been supported in part under DOE grant DE-FG02-91ER40676. 



References 

[1] D. B. Kaplan, Phys. Lett. B288 (1992) 342. 
Y. Shamir, Nucl. Phys. B406 (1993) 90. 



T. Blum et al., |hep-lat/0007038| . 



[2] R. Narayanan and H. Neuberger, Nucl. Phys. B443 (1995) 305. 

[3] H. Neuberger, Phys. Lett. B417 (1998) 141. 
H. Neuberger, Phys. Lett. B427 (1998) 353. 

[4] H. Neuberger, Phys. Rev. Lett. 81 (1998) 4060. 

R. G. Edwards, U. M. Heller and R. Narayanan, Nucl. Phys. B540 (1999) 457. 
A. Borici, Phys. Lett. B453 (1999) 46. 

S. J. Dong, F. X. Lee, K. F. Liu and J. B. Zhang, Phys. Rev. Lett. 85 (2000) 5051. 
T. DeGrand, Phys. Rev. D63 (2001) 034503. 



27 



[5] J. Kiskis, Phys. Rev. D15 (1977) 2329. 

N. K. Nielsen and B. Schroer, Nucl. Phys. B127 (1977) 493. 
M. M. Ansourian, Phys. Lett. B70 (1977) 301. 



[6] 
[7] 



[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 

[17 

[18 

[19 

[20 

[21 
[22 

[23 
[24 



J. Schwinger, Phys. Rev. 128 (1962) 2425. 

C. Gattringer and E. Seiler, Annals Phys. 233 (1994) 97. 

S. Coleman, R. Jackiw and L. Susskind, Annals Phys. 93 (1975) 267. 
S. Coleman, Annals Phys. 101 (1976) 239. 

C. Adam, Annals Phys. 259 (1997) 1 and references therein. 

N. D. Mermin, H. Wagner, Phys. Rev. Lett. 17 (1966) 1133; 

P. C. Hohenberg, Phys. Rev. 158 (1967) 383; 

S. Coleman, Commun. Math. Phys. 31 (1973) 259. 

C. Gattringer, I. Hip and C. B. Lang, Phys. Lett. B466 (1999) 287. 

A. V. Smilga, Phys. Rev. D55 (1997) 443. 

P. H. Ginsparg and K. G. Wilson, Phys. Rev. D25 (1982) 2649. 
M. Liischer, Phys. Lett. B428 (1998) 342. 

F. Niedermayer, Nucl. Phys. (Proc.Suppl.) 73 (1999) 105 and references therein. 

T. Chiu, C. Wang and S. V. Zenkin, Phys. Lett. B438 (1998) 321. 

P. Hernandez, K. Jansen and M. Liischer, Nucl. Phys. B552 (1999) 363. 

P. Hasenfratz, Nucl. Phys. B525 (1998) 401. 

S. Capitani and L. Giusti, Phys. Rev. D62 (2000) 114506 and |hep-lat/0011070| . 

R. Narayanan, H. Neuberger and P. Vranas, Phys. Lett. B353 (1995) 507 and 
Nucl. Phys. Proc. Suppl. 47 (1996) 596. 

J. Kiskis and R. Narayanan, Phys. Rev. D62 (2000) 054501. 

F. Farchioni, I. Hip and C.B. Lang, Phys. Lett. B443 (1998) 214. 

F. Farchioni, I. Hip, C.B. Lang and M. Wohlgenannt, Nucl. Phys. B549 (1999) 364. 

S. Chandrasekharan, Phys. Rev. D59 (1999) 094502. 

W. Bietenholz, I. Hip, Nucl. Phys. B570 (2000) 423. 

28 



[25] G. Katz, G. Batrouni, C. Davies, A. Kronfeld, P. Lepage, P. Rossi, B. Svetitsky and 
K. Wilson, Phys. Rev. D37 (1988) 1589. 

[26] L. Giusti, C. Hoelbling and C. Rebbi, Nucl. Phys. (Proc. Suppl.) 83-84, (2000) 896. 

[27] R. Brower, R. Edwards, C. Rebbi and E. Vicari, Nucl. Phys. B366 (1991) 689. 



29 



